Coherence controlled soliton interactions 
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We demonstrate theoretically and subsequently observe in experiment a novel type of soliton 
interaction when a pair of closely spaced spatial optical solitons as a whole is made partially inco- 
herent. We explain how the character of the soliton interaction can be controlled by the total partial 
incoherence, and show a possibility to change the soliton interaction from attractive to repulsive, or 
vice versa, near a certain threshold in the coherence parameter. 
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In recent years, the study of solitons has attracted con- 
siderable interest in nonlinear optics as well as in other 
areas of physics such as Bose-Einstein condensates [lj. 
Usually, most features of solitons are associated with 
their coherent nature, and the soliton parameters such 
as the amplitude, phase, and frequency are well defined, 
being also responsible for the strength and character of 
the soliton interactions in nonlinear media. 

However, the fundamental concept of solitons as fully 
coherent objects was extended in 1996 to cover more gen- 
eral classes of self-trapped beams which are partially co- 
herent, and therefore can be termed as incoherent soli- 
tons [3 • Somewhat related concepts have been developed 
in the theory of self-trapping of quasi-periodic waves Q 
and spontaneously generated temporal solitons in a non- 
linear medium with instantaneous response |4( . 

To support incoherent spatial solitons, the nonlinear- 
ity must be noninstantaneous, such that the medium is 
insensitive to the fast changing phase of the light beam, 
but only can respond to the time-averaged light inten- 
sity. Most experiments on the incoherent spatial solitons 
have been performed in the biased photorefractive ma- 
terial, whose nonlinear response is noninstantaneous due 
to the fact that its nonlinearity comes from the migration 
of space charges |2|. In parallel to the experimental dis- 
covery, several different theoretical approaches |H H, 
(for an overview of different methods, see Ref. [l|) to de- 
scribe incoherent spatial solitons have been developed in 
the same period of time. One of the theoretical meth- 
ods, the coherent density approach 5], is a powerful tool 
for this kind of analysis, and it is used mostly when the 
dynamics of solitons is the primary subject of study. 

Another phenomenon associated with the properties 
of optical solitons that also caught much research atten- 
tion is their particle-like interaction Q. This has been 
demonstrated, in particular, with the bright optical spa- 
tial solitons in self-focusing media. The main features 
and mechanisms of the interaction of solitons as coher- 
ent objects are well known |2|. If two bright solitons are 
mutually coherent, they attract (repel) each other when 
they are in-phase (ir out of phase). With the relative 
phase between the interacting solitons being other than 



zero or 7r, there is energy transfer from one soliton to the 
other, in addition to the repulsive or attractive interac- 
tion. This energy transfer is most significant when the 
soliton relative phase is it/ 2. On the other hand, if the 
interacting solitons are mutually incoherent or both soli- 
tons are partially coherent (the relative phase between 
them varies much faster than the material can respond), 
the soliton interaction is always attractive [Io| . 

In this Letter, we reveal the existence of a new type of 
soliton interaction that is observed when two interacting 
solitons as a whole are made partially incoherent. We 
find that the soliton interaction dynamics can become 
dramatically different from the case when the solitons 
are mutually coherent. The interaction strength can be 
controlled by the total coherence. Most remarkably, the 
interaction may change from attractive to repulsive near 
a certain threshold in the coherence parameter, or vice 
versa. 

In order to shed light on this novel type of the soliton 
interaction and its origin, first we present our numerical 
results, then provide a qualitatively physics intuition to 
explain why this happens. Finally, we demonstrate ex- 
perimentally the coherence controlled soliton interaction 
in a biased self-focusing photorefractive crystal. 

In this paper, we use the coherence density ap- 
proach p and perform numerical simulations of inter- 
acting partially coherent solitons in a planar geometry. 
Propagation of light beams in a slow Kerr- type nonlinear 
medium can be characterized by interaction of many mu- 
tually incoherent components governed by a set of cou- 
pled nonlinear equations, 



.dE n , ^d 2 E n 



0, 



(1) 



where n is the number of component, \E n \ 2 = I is 
the total beam intensity, x and z are the transverse and 
propagation coordinates normalized to their characteris- 
tic values xo and zo, respectively, D = Azo/(4no£o) is 
the diffraction coefficient, A is the vacuum wavelength, 
no is the medium refractive index, and 7 is the effective 
normalized nonlinear coefficient. 



When at the crystal input the light is generated by 
a partially incoherent source, one can associate different 
components with the same coherent field E s (x) , but hav- 
ing different inclination angles nO s : 



E n (z = 0) = y/G{n0 8 )0 8 E 8 {x) exp(inO s xk ), (2) 

where ho = 2txxqUq/\ and S is a discrete step over 
angles (in radians) which should be chosen sufficiently 
small. We consider a Gaussian angular power distribu- 
tion with a width #o, 



G(0) = (v^/v^o) exp (-26 2 /6l) , 



(3) 



so that the beams are fully coherent for 6q — > 0, and 
become less coherent as 0q is increased. 

We study a mutual interaction between solitons which 
are excited by two beams launched in parallel, 

E a (x) = A {sech[(x + d)/W ] + e^sech[(a; - d)/W }} , 

where 2d is the separation between the beams at the in- 
put, Wo defines the beam widths, <p is the relative phase 
between the two beams, and A = ^/2D/jWq is a coher- 
ent soliton amplitude. 

It is known that two equal- amplitude interacting co- 
herent solitons can experience attraction, repulsion, or 
energy exchange, and all such processes are controlled by 
a single parameter, the soliton relative phase (j). As we 
show below, the interaction of partially coherent beams 
depends critically on the coherence parameter 0q. 

In the simulations, which are performed to supple- 
ment the experimental results presented below, we set 
A = 532nm, no = 2.35, Wo = 5/im, 2d = 20/im, and 
the total propagation length of 25 mm. For the in-phase 
interaction (0 = 0), Fig. HJa) shows that the attraction 
between two solitons makes the solitons launched in par- 
allel cross each other periodically when the entire field 
is coherent (0q = 0). We then increase the coherence 
parameter 0q from zero, i.e. make the entire field par- 
tially incoherent. Since this raises the power threshold 
for a soliton formation, we compensate for this effect by 
increasing the strength of the nonlinearity 7. Our sim- 
ulation results show that the crossing period gets longer 
as the field becomes more incoherent. When #0 reaches 
the threshold value of 0.0028, the solitons become prop- 
agating in parallel [Fig. IHd)] and, finally, when #0 is 
larger than the threshold, the solitons apparently repel 
each other [Fig. [He)]. 

For the 7r-out-of-phase solitons {(j) = tt), the situation 
is just opposite. Figures [if f-j) show that the repulsion 
between two solitons can be changed to attraction when 
the incoherence parameter #0 is larger than 0.0022. Also 
for the 7r/2-out-of-phase interaction (0 = 7r/2), after the 
interaction, the situation that the right soliton is a bit 
brighter than the left soliton due to energy transfer when 
the entire field is coherent (0q = 0) can become oppo- 
site when #0 is above a threshold value around 0.0026 
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FIG. 1: Numerical results, (a)-(e) Interaction of in-phase 
solitons (0 = 0) for different values of the coherence parame- 
ter 0o- (f)-(j) Interaction of 7r-out-of-phase solitons (<fi = tt). 
(k)-(o) Interaction of 7r/2-out-of-phase solitons (<fi = tt/2). 



[Fig. njk-o)], though the opposite energy transfer is not 
that obvious. 

As the next step, we provide a theoretical explana- 
tion on how the degree of coherence (#0) can control the 
soliton interaction. We note that due to the symme- 
try properties of Eqs. QJ, the set of components E n can 
be chosen in a number of different ways which provide 
exactly equivalent descriptions [ll[. In particular, it is 
possible to separate out symmetric and anti-symmetric 
contributions by making a unitary transformation 



E n = (E_ n + E n )/V2 = cos(n6 s x)E s (x)^2G(n6 s )6 s 



for n > 0, E = E = E s (x)y/G(O)0 8 , and 
E n = i(E_ n - E n )/V2 = sm(n0 8 x)E 8 (x)y/2G(n0 a )O 8 , 

for n < 0. For a pair of in-phase solitons (</> = 0), the 
components are symmetric, for n > 0, and they are anti- 
symmetric, for n < 0; the situation is reversed for a pair 
of out-of-phase solitons (<j) = tt). This symmetry is pre- 
served during the propagation, and therefore there is no 
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FIG. 2: Sketch of the experimental setup. 



net energy exchange between the in-phase or out-of-phase 
solitons. Thus, the overall effect of attraction or repulsion 
between the solitons depends on a balance between the 
powers of symmetric (P s ) and anti-symmetric (P a s) com- 
ponents. A specific power ratio p = P s /P&s is required to 
observe a stationary propagation of two solitons in par- 
allel, and it can be estimated using an exact two-soliton 
solution p2 of Eqs. (0) as 



po ~ l + 4exp(-2d/W ) 



(4) 



for 2d ^> Wo. At the crystal input, the parameter p 
monotonously approaches the value of +1 from above 
(below) for a pair of in-phase (out-of-phase) solitons as 
the incoherence parameter Oo is increased. However, 
since the partially-coherent input field does not exactly 
match the soliton profile, some energy is radiated away, 
and we found that the parameter p calculated for the spa- 
tially localized modes may slightly increase or decrease 
during the initial propagation stage. Because of this radi- 
ation mechanism, the condition p = po can be satisfied at 
some threshold value of when the solitons propagate 
in parallel. 

We note that the strongest contribution to the cos- 
type soliton components (n > 0) comes from the angles 
around zero, whereas the sin-type components (n < 0) 
are predominantly excited when n0 s c± A/(4d). There- 
fore, the soliton separation 2d should be roughly inversely 
proportional to the threshold value #Q th \ We have per- 
formed additional simulations to verify this prediction 
and found that when the soliton separation 2d is in- 
creased from 20 /im to 24 /im, the threshold value of 
the coherence parameter is indeed reduced from 0.0028 
to 0.0022, for the interaction of in-phase solitons, and 
from 0.0022 to 0.0018, for the 7r-out-of-phase solitons. 

The interaction of two solitons with the initial <j) = tt/2 
phase shift is more complex. In addition to the effects 
of attraction or repulsion, the solitons can exchange en- 
ergy. At the initial propagation stage, the energy flows 
from the right to the left soliton in components E n with 
n > 0, but the flow is opposite for n < because the 
phase difference is effectively changed in the correspond- 
ing components from tt/2 to —tt/2. As the incoherence is 
increased, the combined powers of two types of compo- 
nents become almost equal, and overall energy exchange 
is largely suppressed as observed in simulations. 
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FIG. 3: Experimental results for the interaction of a pair 
of in-phase solitons for (a) coherent light and (b) incoherent 
light, and a pair of 7r/2-out-of-phase solitons for (c) coherent 
light and (d) incoherent light. 



In our experiment , we use the interaction of stripe pho- 
torefractive solitons 13] as the study object. As schemat- 
ically shown in Fig. [2j we pass the laser beam through 
a rotating diffuser to make the light partially incoher- 
ent. The light beam is then split into two parts. We 
then made these two parts to propagate in parallel by 
using a second beam splitter. These two parts are then 
focused one-dimensionally by a cylindrical lens onto the 
front face of a biased photorefractive crystal. The entire 
crystal is illuminated with a background intensity that 
is a bit stronger than the soliton intensity to make the 
nonlinearity close to the Kerr-type nonlinearity but 
not much stronger to avoid the transverse instability ^| . 

We first take out the rotating diffuser to make the en- 
tire field coherent. We also tune the optical path by 
applying a suitable voltage on the piezo-transducer at- 
tached to the back of one turning mirror to make the 
two parts in-phase. We then observe, due to the in-phase 
attraction, these two parts merge into one at the output 
face of the crystal [Fig.Efa)]. We keep the optical paths 
fixed for the entire system, but now insert the rotating 
diffuser to make the total field partially incoherent. We 
also increase the nonlinearity correspondingly, by higher 
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FIG. 4: Experimental results for the interaction of a pair of 
7r-out-of-phase solitons. The soliton separation at the output 
face decreases from the most coherent light, in the case (a), 
to the least coherent light, in the case (d). 

biasing voltage, due to the fact that incoherent solitons 
need higher nonlinearity [2| . As we can see at the output 
face of the crystal [Fig.|3Jb)], these two parts apparently 
repel each other. This indeed confirms our simulation 
that if the entire field is made incoherent enough, the 
attractive interaction can change to a repulsive one. We 
then tune the relative phase to be tt/2. When the entire 
field is coherent, there is strong energy transfer from the 
left to the right solitons [Fig. Etc)]. However, when the 
entire field is partially incoherent enough, we observe the 
left one is a bit brighter than the right one [Fig. Etd)], 
very similar to the simulation results shown in Fig.^Jo). 

Finally, we set the relative phase to be n. We adjust 
the position of the diffuser and the size of the aperture to 
make the field to be as coherent as possible. Figure Ufa) 
shows the repulsion between the two solitons at the out- 
put face of the crystal. We then increase the incoherence 
by opening the aperture at a larger size. As the incoher- 



ence is increased, the separation become smaller and can 
be smaller than the separation at the input face, indicat- 
ing the repulsion have been changed to attraction. 

In conclusion, we have described theoretically and 
demonstrated experimentally a new type of soliton in- 
teraction which is observed when two solitons as a whole 
are made partially incoherent. The interaction strength 
and type, either attraction or repulsion between the soli- 
tons, can be controlled by varying their total coherence. 
Even more, the interaction may change from attractive 
to repulsive near a certain threshold in the coherence 
parameter. We believe this novel type of the coherence 
controlled soliton interaction is generic for the interac- 
tion of partially coherent waves in nonlinear media, and 
it can be found in other fields. 
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